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We describe a low-temperature approach to the two-dimensional half-filled Hubbard model which 
allows us to study both antiferromagnetism and single-particle properties. This approach ignores 
amplitude fluctuations of the antiferromagnetic ( AF) order parameter and is valid below a crossover 
temperature Tx which marks the onset of AF short-range order. Directional fluctuations (spin 
waves) are described by a non-linear sigma model (NLuM) that we derive from the Hubbard model. 
The parameters of the NLcrM, the spin stiffness and the spin- wave velocity, are calculated as a 
function of the Coulomb repulsion U. The NLctM is solved by a saddle-point approximation within 
the CP^ representation where the Neel fleld is parametrized by two Schwinger bosons. At zero 
temperature, there is always Bose condensation of the Schwinger bosons, which signals AF long- 
range order for any value of the Coulomb repulsion. At finite temperature, the AF long-range order is 
suppressed (in agreement with the Mermin- Wagner theorem), but the AF correlation length remains 
exponentially large. In the CP^ representation, the fermion field is naturally expressed as the product 
of a Schwinger boson and a pseudo-fermion whose spin is quantized along the (fluctuating) Neel field. 
This allows us to write the fermion Green's function as the product (in direct space) of the Schwinger 
boson propagator (which is derived from the NLcrM) and the pseudo-fermion propagator. At zero 
temperature and weak coupling, our results are typical of a Slater antiferromagnet. The AF gap is 
exponentially small; there are well-defined Bogoliubov quasi-particles (QP's) (carrying most of the 
spectral weight) coexisting with a high-energy incoherent excitation background. As U increases, the 
Slater antiferromagnet progressively becomes a Mott-Heisenberg antiferromagnet. The Bogoliubov 
bands evolve into Mott-Hubbard bands separated by a large AF gap. A significant fraction of spectral 
weight is transferred from the Bogoliubov QP's to incoherent excitations. At finite temperature, 
there is a metal-insulator transition between a pseudogap phase at weak coupling and a Mott- 
Hubbard insulator at strong coupling. Finally, we point out that our results straightforwardly 
translate to the half-filled attractive Hubbard model, where the q = (tt, tt) charge and q = pairing 
fiuctuations combine to form an order parameter with SO (3) symmetry. 

PACS numbers: 71.10.Fd,71.10.Hf,71.27.+a 



I. INTRODUCTION 

The Hubbard modeliiSi^ and its generalizations play a 
key role in the description of strongly correlated fermion 
systems such as high-Tc superconductors, heavy fermions 
systems, or organic conductorsi^ Despite its simplicity 
(the model is defined by two parameters, the inter-site 
hopping amplitude t and the local Coulomb interaction 
U, and the symmetry of the lattice), exact solutions or 
well-controlled approximations exist only in a few special 
cases like in one-dimension^ (ID) or in the limit of infinite 
dimension^ 

It is now well established that the ground state of the 
half-filled Hubbard model on a cubic or square lattice has 
antiferromagnetic (AF) long-range orderp2*2i In the weak- 
coupHng Hmit {U <^ it), a Fermi surface instability gives 
rise to a spin-density- wave ground state as first suggested 
by SlaterfiS The AF long-range order produces a gap 
in the quasi-particle (QP) excitation spectrum so that 
the system becomes insulating below the AF transition 
temperature. In the strong-coupling regime {U ^ 4i), 
fermions are localized by the strong Coulomb repulsion 
(Mott-Hubbard locaHzation), thus creating local (mag- 
netic) moments on the lattice sites that are well described 



by the Heisenberg modeli^iii These local moments order 
at low temperature and give rise to a Mott-Heisenberg 
antiferromagnet . 

The main difference between Slater and Mott- 
Heisenberg antiferromagnets lies in the existence or ab- 
sence of preformed local (magnetic) moments above the 
Neel temperature Tn4i In the weak-coupling limit, we 
expect a Fermi liquid phase down to temperatures very 
close to Tn where critical AF fiuctuations start to grow. 
In the strong-coupling Hmit, the system is insulating 
both above (Mott-Hubbard insulator) and below (Mott- 
Heisenberg antiferromagnet) the Neel temperature. 

This simple view, while correct in 3D, breaks down in 
2D. In 2D systems, thermal (classical) fiuctuations pre- 
clude a finite-temperature AF phase transition, and the 
phase transition occurs at Tn = in agreement with 
the Mermin- Wagner theoremii^ Nevertheless, below a 
crossover temperature Tx , the system enters a renormal- 
ized classical regime where AF fiuctuations start to grow 
exponentially. Below Tx, the Fermi-liquid description 
breaks down even at weak coupling, although the sys- 
tem remains metallic. Instead of well-defined Landau's 
QP's, the fermion spectral function A{'k,uj) exhibits two 
(broadened) peaks separated by a pseudogap. 
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The existence of a pseudogap at weak coupling is best 
understood by considering the zero-temperature Hmit. 
At zero temperature, .4(k, uj) is expected to exhibit two 
peaks corresponding to the BogoUubov QP's as in the 
Hartree-Fock (HF) theory. These two peaks are sepa- 
rated by the AF gap which is due to the presence of 
magnetic long-range order. At any finite temperature, 
the AF long-range order disappears in 2D. However, by 
continuity, the two-peak structure in -4(k, u) cannot dis- 
appear as soon as we raise the temperature. As pointed 
out in Ref. 0, the only possible scenario is that at fi- 
nite by low temperature, the fermion spectral function 
exhibits two broadened peaks, which are precursors of 
the zero-temperature Bogoliubov QP's, separated by a 
pseudogap. At strong coupling, the zero-temperature gap 
survives at finite temperature since the system is a Mott- 
Hubbard insulator. 

The simplest description of the AF ground state of 
the 2D half-filled Hubbard model is based on the HF 
theory. It is known that the HF theory remains mean- 
ingful even at large U. In particular, spin-wave modes 
obtained from the Heisenberg model with an exchange 
coupling J — At'^/U can be reproduced from a random- 
phase-approximation (RPA) calculation about the AF 
HF solution li^iiSiiSiil The infiuence of the spin-wave 
modes on the fermionic excitations has been studied 
within one-loopi^ and self-consistent one-loopiSiSfi ap- 
proximations. A QP picture for the coherent motion of a 
particle or a hole appears to be still valid. However, AF 
quantum fiuctuations lead to a significant reduction of 
the Bogoliubov QP spectral weight, with a concomitant 
redistribution of spectral intensity into incoherent excita- 
tions, and a strong renormalization of the AF gap. These 
conclusions are supported by numerical work on the Hub- 
bard model2ii22i2^ and, in the strong-coupling limit, by 
analytical or numerical analysis of the t-J modeliS4 

In spite of its success at zero temperature, the HF the- 
ory fails in 2D since it predicts AF long-range order at 
finite temperature. In the weak-coupling limit, alterna- 
tive approaches, which do satisfy the Mermin- Wagner 
theorem, have been proposed: Moriya's self-consistent- 
renormalized theoryji^i2Si2& the fiuctuation exchange ap- 
proximation (FLEX) (21 or the two-particle self-consistent 
theory»ii None of these approaches gives a unified de- 
scription of the magnetic properties of the 2D Hubbard 
model at finite temperature, both at weak and strong 
coupling. At strong coupling, in the Mott-Hubbard insu- 
lating state, spin degrees of freedom are usually described 
by the Heisenberg model for which various methods are 
available»i2a 

Beside their limitation to the weak-coupling regime, 
these approaches are also unable to account for the 
strong suppression of the amplitude fiuctuations of the 
AF order parameter at low temperature and therefore 
essentially describe Gaussian spin fiuctuations. Be- 
low the crossover temperature Tx, amplitude fiuctu- 
ations are indeed frozen and only directional fiuctua- 
tions [i.e. (transverse) spin waves] survive at low en- 



ergy. The calculation of the single-particle Green's 
function usually relies on a paramagnon-like self-energy 
describing free fermions that couple to Gaussian or- 
der parameter fiuctuationsii^iSL^SiSiiSS This kind of ap- 
proach has been originally introduced by Lee, Rice 
and Anderson to explain the suppression of the den- 
sity of states associated with order parameter fiuc- 
tuations near a charge-density-wave instabilityi^ It 
has been since studied by many authors, in one and 

two dimensi0nS.^^i-^'^i^1^.-^7.^f^r^Pl4(],41,42,4^,44,4R,4f^,47 Tj^g 

sumption of Gaussian spin fiuctuations leads to an overes- 
timation of the fermion density of states at low energy4^ 
Moreover, the artificial presence of amplitude fiuctua- 
tions does not allow to reach the correct T ^ limit i^fi^ 
The effect of "directional" (i.e. phase) fiuctuations of a 
complex order parameter on the fermion density of states 
has been studied both for incommensurate ID Peierls 
system a'^fti'^^ and 2D superconductorsj^ Transverse spin- 
wave fiuctuations in the finite-temperature 2D Hubbard 
model have not received as much attention so far. 

On the experimental side, antiferromagnetism and 
pseudogaps are ubiquitous in low-dimensional strongly 
correlated fermion systems. Pseudogaps were first 
observed in quasi-lD systems near a charge-density- 
wave instabilityi^2i5i More recently, a pseudogap has 
been observed in the metallic phase of high-T^ 
superconductorsi^^i^ Whether the pseudogap in these 
systems is of magnetic or pairing origin is still a mat- 
ter of intense debate. 

In this paper, we describe a theoretical approach which 
provides a unified view of the 2D half-filled Hubbard 
model at low temperature (including T = 0) and for 
any value of the Coulomb repulsioni^ It is based on a 
non-linear sigma model (NLcrM) description of spin fiuc- 
tuations. At zero temperature, our theory describes the 
evolution from a Slater {U ^ At) to a Mott-Heisenberg 
{U :s> At) antiferromagnet. At finite temperature, it pre- 
dicts a pseudogap at weak-coupling due to strong AF 
fiuctuations, and a Mott-Hubbard gap at strong coupling. 
Since it takes into account only directional fiuctuations of 
the AF order parameter, it is valid for T <C Tx, where Tx 
is a crossover temperature which marks the onset of AF 
short-range order. In Ref. one of the present authors 
reported a calculation of the fermion spectral function 
in the weak-coupling limit of the Hubbard model using 
a NLcrM description of spin fiuctuations. However, the 
limitations encountered by previous approaches could not 
be overcome. 

As first shown by Schulz^fi^. spin fiuctuations in the 
2D Hubbard model at low temperature can be described 
by a NLcrM for any value of the Coulomb repulsioni^ 
In Sec. m we give a detailed derivation of the NLcrM 
starting from the Hubbard model. The parameters of 
the NLcrM, the bare spin stiffness and the spin- wave 
velocity c, are calculated as a function of the ratio U/t. 
For U ^ At, we recover the NLcrM derived from the 
Heisenberg model with an exchange coupling J — At'^/U. 
In Sec. mil we introduce the CP^ representation of the 
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NLcrM where the Neel field (giving the direction of the 
local AF order) is expressed in terms of two Schwinger 
bosons. This allows a simple saddle-point solutionii from 
which we obtain the magnetic phase diagram of the 2D 
Hubbard model. At zero temperature, there is conden- 
sation of the Schwinger bosons for any value of U, which 
signals the presence of AF long-range order. At finite 
temperature, the system is disordered by thermal fiuctu- 
ations, but the AF correlation length remains exponen- 
tially large below a crossover temperature Tx (renormal- 
ized classical regim^). In Sec. lIVL we study the fermion 
spectral properties. The fermion is written as the prod- 
uct of a Schwinger boson and a pseudo-fermion whose 
spin is quantized along the (fiuctuating) Neel field. Such 
a decomposition is reminiscent of slave-bosonSS. or slave- 
fermion^*'"^^"®^ theories.-- It allows us to approximate the 
fermion Green's function by the product (in direct space) 
of the Schwinger boson propagator (which is obtained 
from the NLcrM) and the HF fermionic propagator. At 
weak-coupling (U <C 4t) and zero temperature, our re- 
sults clearly describe a Slater antiferromagnet. The AF 

gap 2Ao ~ ie~^'^V^*/^ is exponentially small. As in the 
HF theory, there are well-defined BogoHubov QP's. How- 
ever, because of AF quantum fluctuations, their spec- 
tral weight is reduced by a factor no (0 < no < 1) 
which is given by the fraction of condensed Schwinger 
bosons in the ground state. The missing weight (1 — no) 
is transferred to incoherent excitations at higher energy 
(1 - no < 1 when [/ < 4i). As U increases, the AF 
gap increases and spectral weight is progressively trans- 
ferred from the Bogoliubov QP's to the incoherent ex- 
citation background. At strong coupling {U ^ 4i), our 
results are typical of a Mott-Heisenberg antiferromagnet. 
The AF gap 2Ao is of order U. The incoherent excita- 
tion background carries a signiflcant fraction of spectral 
weight (i.e. no and 1 — no are of the same order) and 
extends over an energy scale of order J — At'^/U above 
the Bogoliubov QP energy ±i?k. At finite temperature, 
the Bogoliubov QP's disappear (no = in the absence of 
Bose condensation) and only incoherent excitations sur- 
vive. Nevertheless, precursors of the zero-temperature 
BogoHubov QP's show up as sharp peaks at ±_Ek in the 
spectral function A(k,uj), with a width of order T. We 
show that these peaks continuously evolve into the zero- 
temperature BogoHubov QP peaks as T ^ 0. This en- 
sures that the spectral function .4(k, u) is continuous at 
the Tn = phase transition. The high-energy incoherent 
excitation background is little affected by a finite temper- 
ature, but the presence of thermal AF fiuctuations gives 
rise to fermionic states below the zero-temperature AF 
gap Ao. At weak coupling, the gap is completely filled 
and replaced by a pseudogap. At strong coupHng, the 
zero-temperature gap survives at finite temperature and 
the system is a Mott-Hubbard insulator. 

On the basis of a numerical calculation in the frame- 
work of the dynamical cluster approximation, Mouk- 
ouri and Jarrell have called into question the exis- 
tence of a Slater scenario in the 2D half-filled Hubbard 



modeliSa*SSiS£ They argue that the system is always a 
Mott-Hubbard insulator at low (but finite) temperature 
even at weak coupHng. We will show that their results 
are not in contradiction with a Slater scenario at weak 
coupling, but merely refiect the strong suppression of the 
density of states due to the pseudogap fSec. lIVCl . 

At half-filHng, the repulsive Hubbard model can be 
mapped exactly onto the attractive model by a canoni- 
cal transformationiS This transformation maps the q — 
(tt, tt) spin correlations of the repulsive model onto the 
q = pairing and q = (tt, tt) charge correlations of the 
attractive model, but leaves the single-particle Green's 
function and the spectral function ^(k, w) invariant. 
Thus the results obtained in this paper apply also to the 
attractive Hubbard model, but with a different physical 
meaning (Sec. EJ. At zero temperature, there is super- 
conducting and charge-density-wave long-range orders. 
As the attractive interaction strength increases, there is a 
smooth crossover from a BCS to a Bose-Einstein behav- 
ior. At finite temperature, the weak-coupling pseudogap 
is due to strong pairing and charge fiuctuations, whereas 
the strong-coupHng gap is a consequence of the presence 
of preformed particle-particle pairs. 



II. DERIVATION OF THE NLaM 

The Hubbard model is defined by the Hamiltonian 

H = -^cl^{i+ fl)Cra + U X^^rT'^rTCr^CrX, (1) 
r.cr r 

where t is the nearest-neighbor hopping operator: 

tCra — t {cic^Q^a -|- Cr— e^cr ^r+eyO" ^ Cr_eyCr) ■ (2) 

At half-filling the chemical potential /i equals U/2. Cx 
and Cy denote unit vectors along the x and y directions. 

(crcr) creates (annihilates) a fermion of spin a at the 
lattice site r. We take the lattice spacing equal to unity 
and set fi = fee = 1 throughout the paper. 

We can represent the partition function of the system 
as a path integral over Grassmann fields Vro-, i'ra- The 
action can be written as S'kin + •S'lnt with 

5kin = / drV^'];(5, -Ai-i>r, (3) 

Sint = U dry^lpr-^A^i^rlAl, (4) 

Jo 

where /3 = l/T is the inverse temperature. In the ki- 
netic action S'kin we have used the spinor representa- 
tion 5" = (V't, V'i)"^- To describe collective spin and 
charge fiuctuations, we introduce auxiliary fields. The 
standard approach is to write the interaction part of 

the action as ^Pt^A^ipt^^ri = j (^'J^'r)^ - j (^'^c^3*r)^ 
and to perform a Hubbard-Stratonovich transformation 
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by means of two real auxiliary fields Acr and Agr- Al- 
though this procedure recovers the standard mean-field 
(or HF) theory of the Neel state within a saddle-point 
approximation, it leads to a loss of spin rotation invari- 
ance and does not allow to obtain the spin-wave Gold- 
stone modes. Fluctuations of and A^r correspond 
to gapped amplitude modes. Alternatively, one could 
write ^int in an explicitly spin-rotation invariant form, 

e.g. VrtV'rTV^ri^ri = -^{"fltT^rf (t - (di, (72,^3) 

denotes the Pauli matrices), and use a vector Hubbard- 
Stratonovich field. Such decompositions, however, do not 
reproduce the HF results at the saddle-point level. ®^ As 
noted earlierjSiS^ this difficulty can be circumvented by 
using the decomposition 

rr^Mn^^rl - \ {^l^rf ^ J {^l^T ■ f2r*r)' , (5) 

where Hr is a site- and time-dependent unitary vector. 
Spin-rotation invariance is made explicit by performing 
an angular integration over fir at each site and time 
(with a measure normalized to unity). The Hubbard- 
Stratonovich transformation then reads 

Eq. Q corresponds to an "amplitude-direction" repre- 
sentation, where the magnetic order parameter field is 
given by A^rJ^r. The HF theory is now recovered from a 
saddle-point approximation over the auxiliary fields Acr, 
Asr and fir (Sec. Ill A|l . Spin- wave excitations can then 
be obtained by considering small fiuctuations of the fir 
field about its saddle-point value. In Sec. IHII we show 
that the amplitude-direction representation ^ allows to 
go beyond the Neel-ordered HF state and derive an ef- 
fective action for the fir field. 



A. HF theory 

Making the ansatz of an antiferromagnetic order, i.e. 
taking constant auxiliary fields A^ = Ac, A^r = A and 
a staggering vector fir — {—lyuz parallel to the z axis, 
one obtains the saddle-point equations 



*Ac = ^^i^l^r), 

A = |(-l)"-(*]:^3*r) 



(7) 

(8) 



At half-filling, the saddle-point value zAc = — C//2 can- 
cels the chemical potential term in ^ . The HF action is 
quadratic, 

Shf= rfrVvft (5, -£-(-!)■- Af73)*r, (9) 

•^0 



and can easily be diagonalized. Due to the translational 
symmetry breaking there is unit cell doubling. In the 
reduced Brillouin zone scheme (|fca;| -I- |/cy| < tt) ele- 
mentary excitations are exhausted by t wo bands of Bo- 
goliubov QP's at energies ±Ek — ±y/ A"^ + ej^, et = 
—2t{cosk X ~\~ cos ky) being the energy of free fermions. 

Using the HF action (pj one obtains the HF single- 
particle Green's function 



^?."^(k,k',c.) 



'5.,c."5,,.'^"^(k,k',c^). 

— Ok.k' , , ^9 + Ok.k' 



(10) 



crA 



(11) 



where tt ^ (tt, tt) and u; = (2n + 1)ttT {n integer) is 
a fermionic Matsubara frequency. The propagator lfT?Hl 
makes it in turn possible to give an explicit form to the 
gap equation JSJ: 



1 

U 



tanh(^) 
2K 



(12) 



We use the notation = J^^ J^. ff^- Eq. JEJ pre- 
diets a phase transition at a finite temperature , 
which is exponentially small at weak coupling and ap- 
proaches U/i at strong coupHng. Similarly to the tran- 
sition temperature, the zero-temperature gap Aq tends 
to U /2 at strong coupling and is exponentially small at 
weak coupling: 
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B. Spin fluctuations 



(13) 



In 2D, the HF theory breaks down at finite tem- 
perature, since it predicts AF long-range order below 
T^^ . Nevertheless, the HF transition temperature bears 
a physical meaning as a crossover temperature below 
which the amplitude of the AF order parameter takes a 
well-defined value. This is sometimes interpreted as the 
appearance of local moments with an amplitude Aq/U. 
Note that at weak coupling the "local" moments can be 
defined only at length-scale of order ~ t/^o, which 
corresponds to the size of bound particle-hole pairs in 
the HF ground state. Thus, stricto sensu, local moments 
form only in the strong-coupling limit when '^-^ 1- 

Below T^^, it is natural to neglect amplitude fiuctua- 
tions of the AF order parameter and derive an effective 
action of the fir field by integrating out the fermionic 
degrees of freedom. We call Tx the crossover temper- 
ature below which AF short-range order appears. As 
will be shown subsequently, in the weak coupling limit 
Tx ~ , whereas at strong coupling Tx ^ J — 
At^/U < T^^. For T < Tx, the amplitude of the 
AF order parameter can be approximated by its zero- 
temperature HF value Aq. Following HaldanejiiiS^ in 
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the presence of AF short-range order {T <Tx) we write 



= (-l)'-nrv/r^L2+Lr 



(14) 



n is a unitary vector representing the Neel field, whereas 
L is the canting vector, orthogonal to n, taking account 
of local ferromagnetic fluctuations, n is assumed to be 
slowly-varying and L to be small. We perform at each 
site and time a rotation in spin space and introduce a 
new fermionic fleld <i>r deflned by — Rr^r- Rr is a 
time- and site-dependent SU(2)/U(1) matrix satisfying 



RrasRl 



(15) 



The above definition means that TZr, the SO (3) element 
associated to Rr, maps onto rir. The U(l) gauge 
freedom is due to rotations around the z axis, which do 
not change the physical state of the system. The spin of 
the pseudo-fermions <I>r is quantized along the rir axis. 
In order to express the action in terms of the new spinor 
variable, it is convenient to make use of the SU(2) gauge 
field Afj^r ~ J2iy=i 2 3 ^/Ir'^i' defined as 



^Or — —RldrRr, 

AfiT = iR\d^Rr, 



^i = x,y. 



(16) 
(17) 



We also define the rotated canting field Ir = T^-^^Lr. 
Given that T^-^^rir = and Lr -L rir, the Ir vector Hes 
in the x — y plane. Using the identity 



(18) 



we re-express the kinetic and interaction parts of the ac- 
tion as 

r/3 



5kin = / dTV$t[a,-Aor 
"'0 



-2t V C0s(-i9^ - A^r)]$r 



(19) 



^int = -Ao rdrV<i>t[(-l)V3v/r^ 
Jo 

+ lr -Crl^r. (20) 

In the above expressions, both 1 and are small, 
since the gauge field is of the order of d^n. We expand 
Eqs. Ijl9l2flll to second order in these variables. To ze- 
roth order, we recover the HF action ^HFi^"] defined in 
(jnj. The first- and second-order corrections in A^^ yield 
paramagnetic and diamagnetic terms S-p and Sd, respec- 
tively. The corrections in I give first- and second-order 
ferromagnetic fluctuations Si and Sp liSi 

Sp ^ - TdT J2 J^r^Mr- (21) 

1^ = 1,2,3 

t r'^ 

= 2 / ^ A;;r'*^cos(-^a^)$r + c.c.(22) 



Si = -Ao 





„=1.2 



Si. = ^ lyrY^i^mij, 



JOt 



The spin-density currents are defined by 



(23) 



(24) 



(25) 



f = <*isin(-ia^)cr^$r H-cc, ^ = x,y. (26) 



We now derive an effective action for the spin variables 
n and L by integrating out the fermions. Keeping terms 
up to second order in A'^^ and 1, the effective action is 
given by first- and second-order cumulants of the four 
perturbative terms S'p, Sd, Si, S12 with respect to the 
HF action: 

^eff[n,L] = {Sp) + {Sd) + {Si) + {S12) 

-'^{Sl),-'^{Sf\-{S,Si)^. (27) 

Evaluation of the first-order cumulants is straightfor- 
ward. Defining Cc as the absolute value of the (negative) 
kinetic energy per site in the HF ground state we have 



(5d) = 



2Ao 
U 



Jo 



2 



{Si) = 0, 



{Si^) - dr 



(28) 

(29) 

(30) 
(31) 



We recognize in equation l|28|l the usual Berry phase 
term. Since it is beheved to play no role in a two- 
dimensional antiferromagneipii we will ignore it in the 
following. 

The calculation of the second-order cumulants seems 
cumbersome at first sight, since it involves (after mov- 
ing to the Fourier space) the current-current correlation 
function H^;^' (q, w^; q', w^,) = (j;;(q, w^)j;;,'(q', w[,))hf- 
[uji, = v2'kT (v integer) is a bosonic Matsubara fre- 
quency.] In fact, as the correlator stands in front of 
second-order quantities, we are interested only in its zero- 
frequency, zero-momentum value IIJ'j^, (0, 0; 0, 0) which 

we denote by HJ^^, . With the exception of HjJ = IIqo 
and nj^^ = Ily^ all these quantities vanish (see Appendix 
El; so that we obtain 

{SD, = f'drY^A^,^' 

JO 

+ r dr J2 Al^\ (32) 
Jo 
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(Sf)^ ^ AlUll r drY^lf, (33) 

JO 

iy = l,2 

{SpSi}^ = AoUll f dr J2 AoX- (34) 
Jo 



Using the invariance of the current-current correlation 
function with respect to rotations of axis , one can es- 
tabUsh the identity ec/2 — nj^^ = 0. It impHes that the 
Af,^ terms in the first- and second-order cumulants can- 
eel each other, which ensures the U(l) gauge invariance. 
The only remaining correlator is the transverse spin sus- 
ceptibility IIoJ = x^. In order to express the effective 
action S'cff[n, L] [Eq. (j27ll ] in terms of n and L we use 
the relations (see Appendix 



1^=1,2 
1.2 



(35) 



(36) 



with C = 1 ioT fj, = x,y and C = — 1 for = 0. Putting 
everything together, we obtain the effective action 



ScS 



1 r^^ 



li=x,y 



+A2(| - x^)l2 - zAox^(nr A hr) • LrjsT) 

where n = drn. Integrating out the canting field with 
the constraint Lr -L rip, we eventually obtain a NLcrM 
for the Neel field: 

^NL.M[n] = ^ r dr I d\ {d^n^y 

(38) 

where we have taken the continuum limit in real space. 
The bare spin stiffness p° and the spin-wave velocity c 
are given by 



1 



(39) 



Eq. l(38ll must be supplemented with a cutoff A in mo- 
mentum space. In the strong-coupling limit, where AF 
fiuctuations are due to local moments, the cutoff A can 
be taken of the order of unity. In the weak-coupling Hmit, 
the Neel field is ill-defined at length-scales smaller than 
^0 ^ t/Ao, since below "local" moments cannot be 
defined (see the discussion at the beginning of Sec. lIIB(l . 
We therefore choose A ~ min(l, 2Ao/c)j2i 

For numerical computation of the spin-wave velocity 
and spin stiffness we use the expressions 



X 



2Al 



4t 



= 7 



^■pa{e)de 



VAT 



(40) 
(41) 



c/t 




FIG. 1: Spin- wave velocity c (solid line) and bare spin stiffness 
pi (dot-dashed line) vs U. For f/ 3> 4f we recover the results 
obtained from the Heisenberg model with J — At^ /U (dashed 
line: spin- wave velocity, dotted line: bare spin stiffness). In- 
set: fraction of condensed bosons at T = 0. 1 — no is expo- 
nentially small at weak coupling {U <^ At), while no — 0.6 
for U ^ At. no determines the mean-value of the Neel field 
in the ground state ((nr) = uqUz) and the spectral weight of 
the Bogoliubov QP's (see Sec. IIV|l . 



po is the density of states of free fermions on a square 
lattice. It can be expressed, using the complete elliptic 
integral of the first kind K, as po{e) = {2T:'^t)~^K[{l — 
e^/16i^)^/^] for |e| < At. In the strong coupling limit we 
recover the results obtained from the Heisenberg model 
with an exchange coupling J — At'^/U: the spin stiffness 
equals J/4 and the spin-wave velocity \/2J. At weak 
coupling, c goes to zero like 2T:^^^'^t{U/t)i and pi*^' ~ t. 
The factor {U/t)^ is due to the Van-Hove singularities. 
The results are shown in Fig.^ 



The NLcrM defined by Eqs. II38I41|I was first obtained 
by SchulziS The value of the spin- wave velocity agrees 
with the result obtained from an RPA calculation about 
the zero-temperature AF HF statei ^^i^^i^'^ 



Note that we also expect a damping term with a char- 
acteristic frequency Ugf in the NLcrM action (j38ll at weak 
coupling. This term comes from the damping of spin 
fiuctuations by fermion excitations which are gapless in 
the weak coupling Hmit (see Sec. lIVII . It is missed in our 
approach since we expand around the zero-temperature 
AF state which has only gapped quasi-particle excita- 
tions. In the renormalized classical regime, this term is 
however negligible since cugf « (critical slowing 

down),i^ 
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III. MAGNETIC PHASE DIAGRAM 



T 



Let us recast the NLcrM action in a more usual form, 
by making use of the couphng constant g ~ c/p'^: 



'5'NLcrM[n] — 



1 

2^9 



/3 



dr / d^r 



1 



ti=x,y 



(42) 

To solve the NLtrM, we use a saddle-point approximation 
in the CP^ representation^ii which proves well suited for 
the computation of the fermion Green's function. In the 
CP^ representation, the Neel field is expressed in terms 
of two Schwinger bosons. 



rir = zlcrz 



(43) 



with Zr = (zrTj^ri) • The condition nf — 1 translates 
into zjzr = 1. The rotation matrix R can be expressed 
as 



Rr = 



-Z, 



(44) 



The U(l) gauge symmetry now manifests itself in the 
invariance of the rir vector and the relation ltT5|l defin- 
ing the rotation matrix under the transformation Zr 
e*"--Zr. 

The NLcrM expressed in terms of Schwinger bosons 
involves terms quadratic and quartic in Zr. The latter 
turns out to be proportional to A^r^ with expressed 
in terms of Schwinger bosons [Eqs. Ijl7|l and 14411 ]. It is 
decoupled by an auxiliary field a^r- To handle the uni- 
modularity condition zjzr — 1 one introduces Lagrange 
multipliers Ar at each time and site. The partition func- 
tion then becomes (see Ref. O) 



Z = 
S = 



J V[z,af,,\]t 



(45) 



dr / 



iAr(z]:Zr - 1) H |(9r - aOr)ZrP 

gc 



2c 



fj,=x,y 



(46) 



the Zri and Zri being now unconstrained bosonic fields. 
One then performs a saddle-point approximation over the 
Ar and a^^r fields. When the CP^ representation is gen- 
eralized to the CP^~^ representation by introducing N 
different z bosons, the approximation becomes exact in 
the limit N — > ooiii Within the ansatz of a uniform static 
saddle-point solution iXr — 2m? / gc and a^ir — 0, the 
propagator can be read off from Il4(i|l : 



= 5c^^c^'5uJ,,,uJ'Ja^a'^^a{fi■, Wi.),(47) 



2?<T(q,t^iy 



-9C 



- /37V71o(5<^,T'^a;„,0'5q,0, (48) 

(49) 



QC 



RC 



QD 



- g 



Hubbard model 

FIG. 2: Phase diagram of the NLctM derived from a saddle- 
point approximation in the CP^ representation. At T = 0, 
there is long-range order when the coupling constant g < 
Qc — 47r/A. The three finite-temperature regimes correspond 
to "renormalized classical" (RC), "quantum critical" (QC) 
and "quantum disordered" (QT))-^^ The ground state of the 
2D half-filled Hubbard model is ordered for any value of the 
Coulomb repulsion U. At finite temperature, there are strong 
AF fiuctuations with an exponentially large correlation length 
C > 1 (RC regime). 



where Af is the number of lattice sites. The saddle-point 
equation for the Lagrange multiplier reads 



-T 



q|<A +Wq 



no 



(50) 



In Eqs. Il47l5fl|l . we have allowed for a Bose condensa- 
tion of the Schwinger bosons in the mode q = 0, with 
n-o = (z'f(q = 0, = 0)z(q = 0, = 0)) the frac- 
tion of condensed bosons. Bose condensation signals the 
appearance of AF long-range order: (rir) —noUz. Know- 
ing the propagator of the z field, one can then calculate 
the spin-spin correlation function using Eq. (021 • The 
AF correlation length ^ is related to the mass m of the 
bosonic propagator V via m = c/2£^w^ m vanishes when- 
ever the fraction of the condensed bosons is finite. 

At zero temperature, the solution of the saddle-point 
equation lf5?ill shows that the NLcrM is ordered at small 
g (m = and no > 0) and disordered by quantum fiuc- 
tuations at large g (to > and no = 0) . The two regimes 
are separated by a quantum-critical point at gc = 47r/A. 
In the ordered phase {g < gc) , the fraction of condensed 
bosons is no = 1 — g/gc- 

The condition of zero-temperature long-range order is 
satisfied in the NLcrM derived from the half-fillfid. Hub- 
bard model (Fig.EJ. For U < At, g/gc ~ e-'^^'^Fl^ is ex- 
ponentially small. For U > 4t, p° ~ J/4 and cA ~ %/2J, 
so that g/gc — V2/'k < 1. Notice that setting the cutoff 
to a higher value at strong coupling would lead us into 
the quantum disordered regime. However, our choice is 
consistent with results obtained by mapping the Hubbard 
model at strong coupling onto the Heisenberg model. It 
is known, both from numerical and analytical work, that 
the 2D quantum Heisenberg model on a square lattice is 
ordered at zero temperaturei£^ 



8 



T 









- 






■ 


// . 






',1 R^C regime 

/ : c > 1 


Mott^HuUard y 
insulator 




/ PG : 






huierj. Neel order 





10 20 

U/t 



FIG. 3: Phase diagram of the 2D half-filled Hubbard model. 
T > T^": Fermi liquid (FL) phase; Tx < T < T^^: local 
moments with no AF short-range order (Curie spins, ^ ~ 1); 
T = 0: Slater {U < 4t) and Mott-Heisenberg {U > 4t) 
antiferromagnets. At finite temperature, there is a pseudogap 
phase {U < 4i) and a Mott-Hubbard insulator {U > U) 
separated by a metal-insulator transition (dotted line) defined 
by the vanishing of the tunneling density of states p{ijj = 0) 
at zero energy (Sec. lIVIl . All lines, except Tn = (thick 
solid line), are crossover lines. The NLcrM description is valid 
below Tx (RC regime). [From Ref.Isl.] 

Fig. n] shows the fraction of condensed bosons as a 
function of U. For this, and subsequent, numerical 
calculations we use a smooth cutoff, i.e. i|q|<A ~^ 

/q 'i-e-TOSo"^" • ■'■^ contrast to a hard cutoff, this pro- 
cedure prevents artificial features in the fermion spectral 
function and in the density of states. The parameter 
go is adjusted so as to reproduce in the strong-coupling 
limit {U 3> 4:t) the result | (rir) | = no — 0.6 obtained 
from the Heisenberg modeli^^ While the value of no for 
U ^ At and U ^ At does not depend on ^o, the behavior 
at intermediate coupling is strongly cutoff dependent. 

At finite temperature, the AF long-range order is sup- 
pressed (no = 0, rn > 0), in agreement with the Mermin- 
Wagner theorem. For systems that exhibit AF long-range 
order at T = 0, the correlation length remains neverthe- 
less exponentially large at low temperature [renormalized 
classical regime, see Fig. J^J]. From Eq. H5()ll . we deduce 

C = -^, m = Te-^, (51) 

where ps = Ps(l— g/^c) is the zero-temperature spin stiff- 
ness. The mass m of the bosonic propagator being much 
smaller than the temperature, the dominant fluctuations 
are classical. 

Let us now discuss the Hmits of vaHdity of the NLcrM. 
The derivation of the NLcrM is based on the assumption 
that the dominant low-energy fluctuations are transverse 
spin waves with a large correlation length. The condition 
T ^ T^^ ensures that amplitude fluctuations of the AF 
order parameter are frozen at low energy. One should 



also verify that the computation of ^ within the NLcrM is 
consistent with the assumption of AF short-range order, 
i.e. ^ » A^^ or, equivalently, m ^ cA/2. We define T' 
as the solution of the equation m ~ cA/2 obtained from 
Eq. jnni)- Then, the domain of vaHdity of the NLcrM is 
given by T ^ Tx ^ min(T-^^, T'). At weak coupling, 
Tx ~ ^ while Tx ~ T' ~ J at strong coupling. The 
crossover temperature Tx displayed in Fig. is a smooth 
interpolation between and T'.^'* 

The phase diagram is shown in Fig.El Above , spin 
fluctuations are not important and we expect a Fermi 
liquid behavior. Between T^^ and Tx (a regime which 
exists only in the strong-coupling limit), local moments 
form but with no AF short-range order (Curie spins: 
^ ^ 1). Below Tx, the system enters a renormalized 
classical regime of spin fluctuations where the AF corre- 
lation length becomes exponentially large [Eq. AF 
long-range order sets in at Tn = 0. Although there is a 
smooth evolution of the magnetic properties as a func- 
tion of C/, the physics is quite different for U <^ At and 
U ^ At. This will be shown in Sec. IIVI by studying 
the fermion spectral properties. The main conclusions 
are shown in Fig. El At zero temperature the system 
is an antiferromagnet, which evolves from a Slater to a 
Mott-Heisenberg behavior as U increases. At flnite tem- 
perature there is a pseudogap phase for C/ <C 4t and a 
Mott-Hubbard insulator for U ^ At. These two regimes 
are separated by a (flnite-temperature) metal-insulator 
transition (dotted line in Fig.|S|l deflned by the vanishing 
of the tunneling density of states p{uj = 0) at zero energy. 



IV. FERMION SPECTRAL PROPERTIES 

In this section, we study the influence of the long- 
wavelength spin fluctuations on the fermion spec- 
tral properties. The fermionic Green's function 
— (^riri^J2r2)) written here as a 2 x 2 matrix in spin 
space, can easily be related to the pseudo-fermions by 
use of the relation = Rr^r- 

g(l,2) = -(i?i$ici>ti?t^. (52) 

Here we use the shorthand notations 1 = (ri, n) and 2 = 
(r2, T2)- The averaging in the above expression should be 
performed with respect to the action S'HF[*i'] + 'S'iz, $, L] 
obtained in Sec. lH BI from the second-order expansion in 
L and 9^n. S' stands for the sum of the perturbative 
corrections Sp, Sd, Si, Sp deflned in H21I24|I . Integrating 
flrst the pseudo-fermions, we can write the propagator as 

0(1,2) = 1 y"l?[z]e-^~^'"Wi?i5(l,2|z)i?|, (53) 
Z = / X'[z]e"^"^''"[^l, (54) 
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where 5(l,2|z) is the pseudo-fermion propagator calcu- 
lated for a given configuration of the bosonic field z: 



S(l,2|z) = - 



W2 



*„-Shf[*]-S'[z,*,L] 



(55) 



The action ^hf [^] + <5" [z , $ , L] describes HF fermions in- 
teracting with spin fluctuations via the action S' . Since 
the HF pseudo-fermions are gapped, we expect a per- 
turbative expansion in S" to be well-behaved. To lead- 
ing order, Q{l,2\z) = Q^^{1,2) and the fermion Green's 
function simplifles into 



(56) 



where the product of rotation matrices is averaged with 
the NLcrM action. This approximation neglects the ef- 
fect of spin fluctuations on the propagation of pseudo- 
fermions. Their influence on the propagation of fermions 
is implemented only through the decomposition of the 
fermion into a boson and a pseudo-fermion. 

Using the Schwinger boson propagator derived in 
Sec. Uni [Eqs. ^47149^ ]. we have 

((-Rl)criQ!i (-^2)0-202 ) ^ ^'^(Tl, 0-2^01,02 

X (P(l,2)-5,,,„,no),(57) 

where V is the non-condensed part of V^j. Using this 
expression in Eq. (I56II we flnally obtain for the fermion 
Green's function: 



-(V'rr.VrV.') = 5„ ,a'Q a{r ^ , T - t') , (58) 

e.(k,k',c.) = -%^E / er(k-q,k-q,c.-^.)P(q,c..)+noSHF(k,k',u;). (59) 

Since tiq vanishes at flnite temperature, the fermion Green's function is spin-rotation and translation invariant in 
the absence of AF long-range order. We show below that the flrst term of the rhs of Ij59|l corresponds to incoherent 
excitations. At zero temperature, the last term of l|59ll describes Bogoliubov QP's carrying a total spectral weight tiq. 

To study in detail the fermion excitations, we consider the spectral function ^(k,a;) = — 7r~^ImC/CT(k, k, io; 
Lu + iO+) and the tunneling density of states (DOS) p{uj) = J dujA(k,uj). Performing the summation over bosonic 
Matsubara frequencies in Eq. Ij59|l we obtain 

^(k,cj) ^ Anc(k,t^)-Fno-4HF(k,cj), (60) 

Anc(k,w) = / ;^|[71_B(tJq) + n_F(-£^k-q)] [Uk-q^(^ ^ ^ ^k-q) + I'k-q'^l^ + + ^k-q)] 

+ ["S(t^q) + ^^(i^k-q)] [U^_^5{UJ + Uq - Sk-q) + v'^_^5{uj -UJq + Sk-q)] |, (61) 



where np{Ljj) and nB{(^) are the usual Fermi and Bose oc- 
cupation numbers (e'^" ±1)^^, and -4hf the HF spectral 
function: 



^HF(k, Lu) = uld{Lu - Ek) + vl6{uj + Sk) 



ek 

El, J 



1 



(62) 
(63) 



One can check that the spectral function ^(k, uj) is nor- 
malized to unity. From Eqs. (I6r)l61|l we deduce 



/ dujA{k,uj) = / — (nB{ojq) + ^ 

J J\<i\<A V ^ 



no 



= 1, 



(64) 



where the second equality is obtained by using (zjzr) = 1 
[Eq. (EDI]. From Eqs. H6(163|l . we obtain 



(65) 



Pinc('^) 
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(66) 



Aire 



nB(tj')pHF('^ + uj')9{uj + lu') 



+ (nBiiu') + l)pHFii^ - i^')9{u) - Lo') , (67) 



where 



Phf(c^) =0(^'- A^) 



^0 



Po 



^0 



is the HF DOS and the step function. We have ap- 
proximated the Fermi occupation numbers by their zero- 
temperature limit, which is valid for T <C . 
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FIG. 4: Spectral function -4(k, lo) in the weak-coupling limit 
U — t for T = (Slater antiferromagnet) and T — Ao/5 
(pseudogap phase), k = (7r/2,7r/2). The vertical lines repre- 
sent Dirac peaks of weight no/2 (Bogoliubov QP's). At finite 
temperature, precursors of the zero-temperature Bogoliubov 
QP's show up as peaks of width ~ T at ±i5k. At low energy 
(and T > 0), we observe a pseudogap with an exponentially 
small spectral weight at uo = Q. Energies are measured in 
units of t. [From Ref. I5fli .] 




FIG. 5: Spectral function ■4(k,tj) in the strong-coupling 
regime U = 12t for T = (Mott-Heisenberg antiferromagnet) 
and T = J/5 (Mott-Hubbard insulator). At T = 0, when U 
increases, spectral weight is transferred from the Bogoliubov 
QP peaks to the incoherent excitation background. [Note the 
difference in the energy scale, which is fixed by Aq, between 
Figs.HandlH] [From Ref. 



A. T = 0: Slater vs Mott-Heisenberg 
antiferromagnetism 

At zero temperature, the incoherent part of the spec- 
tral function [Eq. ijfiTll ] can be simplified. All the occupa- 
tion factors vanish, except fermionic factors at negative 
energies which are equal to 1, so that 




+vl_^S{uj + LUq + Ei,^q)]. (69) 
In the same way, we obtain for the DOS 

Pi'nci^) " ^ duj'paF{i^ - uj')0{uj - uj'). (70) 

In Figs. I4I5I we show the spectral function at the 
k — (7r/2,7r/2) point of the non-interacting Fermi sur- 
face at weak (U — t) and strong {U = \2t) coupling. 
The spectral function yl(k, w) exhibits a gap 2Ao, which 
is a consequence of AF long-range order. There are well- 
defined Bogoliubov QP's with excitation energy ±i?k, as 
in HF theory, but their spectral weight is reduced by a 
factor no < 1 because of quantum spin fluctuations. The 
remaining weight (1 — no) is carried by an incoherent 
excitation background at higher energy (|a;| > £^k)- 

There are important differences between the weak 
{U <^ At) and strong ([/ 3> 4i) coupling regimes. First, 

the AF gap 2Ao ~ ie~^^V^*^ is exponentially small at 
weak coupHng, while it tends to C/ for C/ 3> 4t. Second, 
the Bogoliubov QP's carry most of the spectral weight 
in the weak-coupling regime, since g/gc = I — uq is ex- 
ponentially small when ?7 <C 4^. As U increases, spectral 
weight is transferred from the Bogoliubov QP's to the in- 
coherent excitation background, and at strong coupling 



(U ^ At) the incoherent excitation background carries a 
significant fraction of the total spectral weight (i.e. no 
and 1 — no are of the same order). Third, the energy 
range of the incoherent excitation background depends 
on the value of U. From Eq. Ij69|l we see that it extends 
from Ek to \/E^ + IGt'^A'^ + cA. At weak coupling, 
the upper limit turns out to be of order Aq (for k lying 
on the non-interacting Fermi surface). Thus, the energy 
range of the incoherent excitation background remains 
very small with respect to the dispersion of the Bogoli- 
ubov QP energy i?k, which is of order t when Ao <C t. 
At strong coupling, the incoherent excitation background 
above E\s_ ^ U/2 extends over a range of order J. This 
energy range is of the same order of magnitude as the 
dispersion of the Bogoliubov QP energy, which is also of 
order J when C/ >• 4t [as can be seen from the expansion 
E^k ^ U/2 + J(coskx + cos ky)'^]. 

In Figs. I6I7I we compare the zero-temperature DOS 
p(w) and the non- interacting DOS Poi^)- At weak cou- 
pling (U = t), p(tjj) is similar to the HF result, with 
no visible effect of the incoherent excitation background. 
p{uj) differs from pi){Ljj) mainly at low energy, due to the 
(small) AF gap 2Ao. At strong coupling {U — 12i), p(w) 
differs strongly from Pq{uj), due to an AF gap 2Ao ~ U 
exceeding the non-interacting bandwidth. There is also 
a significant difference between p{u}) and pHF(t^), which 
results from the incoherent excitation background. 

The spectral function ^(k, w) and the DOS p{uj) are 
typical of a Slater antiferromagnet at weak coupling and 
of a Mott-Heisenberg antiferromagnet at strong coupling. 
As shown in the next section. Slater and Mott-Heisenberg 
antiferromagnets behave very differently at finite temper- 
ature. 
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FIG. 6: Zero-temperature DOS p{i^) at weak coupling: U — t 
(Slater antiferromagnet). p(<^) differs from the free-fermion 
DOS po(i^) only at low energy due to the opening of the 
AF gap 2Ao (see inset). Since the incoherent excitation 
background carries a negligible fraction of the total spectral 
weight, there is no noticeable difference between p(^^) and the 
HF DOS phf(w) (not shown in the figure). 
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1. Precursors of Bogoliubov QP's 

At finite temperature, the coherent part of the spec- 
tral function disappears. However, sharp peaks are still 
observed at the HF energy ±i?k- To study the peak at 
i?k, let us perform a few approximations on the finite- 
temperature spectral function l(fiT|l . First, at positive 
energies, almost all the spectral weight comes from the 
terms proportional to w^-q in (except at energies 
close to zero), whose sum will be denoted by . Sec- 
ond, we replace the Fermi occupation number by the step 
function, given that the temperature is small compared 
to i?k- Regrouping terms containing the Bose occupation 
numbers we obtain 



A>(k,L0) 

> I 
peak V 



= ^>g(k,C.)+^>,,(k,C^), (71) 
+ 5{U0 + t^q - Sk-q)] . (72) 



.4^^ has the same expression as the incoherent excitation 
background term ijfifljl at zero temperature. It thus de- 
scribes a temperature-independent incoherent excitation 
background at energies above i?k- -^pcak gives rise to the 
peak at the HF energy E\^. To see this, let us put it into 
a more expHcit form. Because of the bosonic occupation 
numbers, the sum over q in Ij72|l is dominated by wave- 
vectors satisfying Wq < T or, equivalently, |q| < T/c. 
For T ^ Tx, T/c <C 1 and we can neglect the q depen- 
dence of i?k-q and Uk-q- The integrand then becomes 
isotropic, and one can use / = /m^"|^- The result 



^poak(k,^) 



2 y /\ 



(73) 



FIG. 7: Same as Fig. but at strong coupling: U — 12t 
(Mott-Heisenberg antiferromagnet). p{iu) differs strongly 
from the non-interacting DOS po(<^), as the AF gap exceeds 
the non-interacting bandwidth. It also differs from the HF 
DOS phf(i^) due to the incoherent excitation background car- 
rying a significant fraction of the total spectral weight. 



B. T > 0: pseudogap vs Mott-Hubbard gap 



At finite temperature, no vanishes and A = Ainc- The 
result of the numerical calculation for U = t and U = \2t 
for k (7r/2,7r/2) is shown in Figs.21and|3 ^(k, w) ex- 
hibits broadened peaks of width T at the HF QP en- 
ergy ±i?k- These peaks are incoherent precursors of 
the zero-temperature Bogoliubov QP peaks. The zero- 
temperature AF gap is partially filled at strong coupling 
and transforms into a pseudogap in the weak-coupling 
regime. At higher energy (|w| > i?k), a roughly feature- 
less incoherent excitation background is observed. 



for \uj — E\i\ > m, and vanishes for \uj — i?k| < m. For 
TO < Iw - £;k| < T, A^^^^{\<i,uj) behaves like T/\uj - 
i?k|- At energies further away from the peak center, it 
decreases like e"!"^"^"*!/^. Thus the width of the peak is 
of the order of the temperature and therefore corresponds 
to incoherent excitations. The vanishing of .4(k, w) for 
1^^ — £^k| < TO is clearly unphysical (note that it cannot 
be seen in the figures, since to is exponentially small). It 
would be suppressed by any finite life time in the bosonic 
propagator V. The finite-temperature DOS suffers from 
the same artifact (i.e. p(a') = for |w — Ao| < m). 
The spectral weight of the peak at i?k is 



C^'^-^pcak 



(k,C.) = ul^T\u(- 

ZTTC \ TO 



9c 



(74) 



where the last result is obtained using Eq. Ij51|l . The 
spectral weight of the peak turns out to be tempera- 
ture independent and equal to u^tiq (uq = 1 — g/gc), 
which is nothing else but the Bogoliubov QP weight in 
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FIG. 8: Finite-temperature spectral function at weak cou- 
pling for two different points of the non-interacting Fermi sur- 
face. For k close to (7r/2,7r/2) a second peak appears below 
Ao (see text). 



the ground state. We conclude that the peak is an inco- 
herent precursor of the zero-temperature Bogoliubov QP 
peak. As the temperature decreases, it retains its spec- 
tral weight but becomes sharper and sharper, and even- 
tually becomes a Dirac peak at T = 0. As expected, the 
spectral function evolves continuously when T ^ 0. As 
in the zero-temperature case, the dependence of no upon 
U describes the transfer of spectral weight from the Bo- 
goliubov QP's to the incoherent excitation background 
when the Coulomb repulsion increases. 

The approximation ifTsjl suggests that the peak in 
^(k, w) should exhibit the same features, regardless of 
the location of k on the non-interacting Fermi surface. 
Numerical calculations confirm this conclusion, with one 
exception. For wave-vectors near (7r/2,7r/2), a second 
(smaller) peak appears at low energy (Fig. IHJ . From a 
mathematical point of view, it is due to the vanishing 
of the first-order derivative of the argument of the delta 
function in Eq. lf72|l . which occurs for VqWq — Vq£'k-q. 
The energy at which the integration contour in the q 
plane, defined by the delta function , passes throu gh this 
point can be estimated to be Ao\/l — (c/|vk|)^, where 
Vk = Vkek is the free-fermion velocity. For wave-vectors 
verifying |vk| < c, i.e. sufficiently close to the Van-Hove 
singularities, the second peak disappears. We believe this 
second peak to be an artifact of our lowest-order approx- 
imation in the pseudo-fermion-boson interaction. 



vs MoU-Hubbard 



As shown in Figs. 14151 the spectral function ^(k, w) 
extends below the HF energy (and above — i?k for lu < 
0) at finite temperature. The corresponding contribution 
to .4(k, uj) is given by [see Eq. (EU] 
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FIG. 9: Finite-temperature DOS p(i^) at weak coupling: U = 
t, T — Ao/5 (pseudogap phase). At lo — the DOS is finite 
but exponentially small [Eq. I|76|l]. 
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FIG. 10: Finite-temperature DOS p{uj) at strong coupling: 
U = 12t, T = J/5 (Mott-Hubbard insulator). 



(75) 



The presence of the Rose occupation number n_B(wq) 
shows that the low-energy fermion states {\uj\ < i?k) are 
due to thermal bosons, i.e. thermally excited spin fiuctu- 
ations. A fermion added to the system with momentum k 
and energy \uj\ < i?k can propagate by absorbing a ther- 
mal boson of energy ujq and emitting a pseudo-fermion 
with energy -Ek-q = uj + ujq. 

The lowest fermion energies are obtained by solving 



(or LU 



-£'k+q + Wq). In the weak 



coupling limit, maxq(Li;q) — cA ^ 2Ao and i?k-q ~ Ei^- 
Thus there is spectral weight at zero energy: the spectral 
function and the density of states exhibit a pseudogap 
(Figs.Elandini. Note that the DOS remains exponentially 
small at low energy: 



cosh 



(?)■ 



\uj\ < Aq. 



(76) 



This result differs from pseudogap theories based on 
Gaussian spin fluctuations which flnd a much weaker 
suppression of the density of states at low energy4^ It 
bears some similarities with the results obtained by Bar- 
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Tjt 




FIG. 11: Lines p(io = 0) = const in the {U,T) plane. The 
vertical line corresponds to p{u: = 0) = 0. 



tosch and Kopietz for fermions coupled to classical phase 
fluctuations in incommensurate Peierls chainsi^ In the 
strong coupHng limit, thermally excited spin fluctuations 
lead to a small reduction of the zero-temperature gap 
since cA J < i?k ^ U/2. The system is a Mott- 
Hubbard insulator with a gap 2Ao of order U (Figs. El 
and [19 . 

A last comment is in order here. Since the system 
is in the renormalized classical regime, it is tempting to 
treat the NLcrM in the classical limit (which amounts 
to neglecting the quantum (temporal) fluctuations of the 
Neel field n) . Such an approach is expected to be at least 
qualitatively correct for the low-energy bosons (wq < T) 
and should then give a good approximation of ^(k, uj) 
in the vicinity of the peaks around u = ±i?k. Retaining 
only the = contribution in Eq. Ij59|l . one finds 



Ai(k,c.)=T / ^ 



(77) 



Eq. Ij77|l can also be obtained from Eq. Ij61|l by using 
nsi^^q) -I- 1 ~ n_B(wq) ~ T/ujq ^ 1 and neglecting the 
term ±Wq in the argument of the delta functions. It 
is readily seen that the classical calculation does not 
reproduce the pseudogap, since ^ci(k, w) vanishes for 
\uj\ < Eu- Although the pseudogap originates from ther- 
mally excited spin fluctuations in the renormalized clas- 
sical regime, a fully quantum-mechanical calculation of 
.4(k, to) turns out to be necessary to account for the pres- 
ence of low-energy fermion excitations. 



C. Finite-temperature metal-insulator transition 



fe rmion to t he maximum energy of a Schwinger boson 
Vto2 -hc2A2. For T ^ the result is Uc ~ 4.25i. It 
should be noted that the NLctM, which is a low-energy 
theory, does not allow us to describe accurately the high- 
energy Schwinger bosons (with |q| ~ A) and in turn the 
low-energy fermion excitations. In particular, the criti- 
cal value of U calculated above depends on the cut-off 
procedure used in the NLcrM. Note also that we do not 
know at which temperature and how the metal-insulator 
transition ends. 

Fig. El shows the fines p{u} = 0) = const in the {U,T) 
plane. Our results are in (semi-quantitative) agreement 
with the numerical calculation of Moukouri et al^ Using 
the criterion p[ijj = 0) < 10-^/(2*) to identify the Mott- 
insulating phase, these authors concluded that the sys- 
tem is always insulating at low (but finite) temperature 
even in the weak-coupling limit, which seems to invali- 
date the Slater scenario as the mechanism for the metal- 
insulator transition (which requires Tmit — = 
Our approach shows that the results of Ref . 01 are not in 
contradiction with a Slater scenario at weak coupfing, but 
merely reflect the exponential suppression of the density 
of states due to the presence of a pseudogap. A similar 
conclusion was reached in Ref. 



V. ATTRACTIVE HUBBARD MODEL 

In this section, we show that the results obtained in 
the previous sections translate directly to the attractive 
Hubbard model. The latter is defined by the Hamiltonian 

H^~J2 ""raii + ^^)Cra ^ U 4^ ^^1; 4^0^,1, (78) 



where —U (C/ > 0) is the on-site attraction, fi = —U/2 
at half-filling. 

Under the particle-hole transformation^ 

Crl - i-iycl^, - (-l)""Cri, (79) 

the Hamiltonian becomes (up to a constant term) 

r.a r 

-(M + C//2)^(ci^CrT-ci^Cri), (80) 

r 

and the charge-density and pairing operators transform 



We conclude from the results of Sec. IIV Bl that our 
approach predicts a finite-temperature metal-insulator 
transition between a pseudogap phase and a Mott- 
Hubbard insulator as the strength of the Coulomb in- 
teraction increases: at a critical value Uc, the density of 
states at zero energy p{uj = 0) vanishes and the pseudo- 
gap becomes a Mott-Hubbard gap (Fig. EJ. Uc is ob- 
tained by equating the minimum energy Aq of a HF 



Pr 



25*;' + 1, 



(81) 
(82) 

(83) 



where Sr = cJ<TCr/2 and = 5"^ ± iS^. The trans- 
formed Hamiltonian Ij8()|l corresponds to the repulsive 
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FIG. 12: Phase diagram of the 2D half-filled Hubbard model 
with an attractive interaction —U {U > 0). T > T^^: Fermi 
liquid (FL) phase; 2x < T < T™: preformed pairs with 
no superfluid or charge-density- wave order ~ 1); T < Tx: 
renormalized classical (RC) regime ^ 1); T = 0: super- 
conducting (SC) and charge-density-wave (CDW) long-range 
orders {U < it: BCS limit, U > At: Bose-Einstein (BE) 
limit). The dotted line is obtained from the vanishing of the 
tunneling DOS p{uj = 0) at zero energy. All lines, except 
Tc = (thick solid line), are crossover lines. 



half-filled Hubbard model with a uniform magnetic field 
/i -f C//2 along the z axis coupled to the fermion spins. 
At half-filling (/i — —U/2), the latter vanishes and the 
Hamiltonian II8()|I reduces to the one studied in the previ- 
ous sections. Thus, in the attractive model, q = tt charge 
and q = pairing fluctuations combine to form an order 
parameter with SO (3) symmetry. Away from half-filling, 
the degeneracy between charge and pairing fiuctuations 
is lifted (by the uniform magnetic field ^ + U/2 in the re- 
pulsive model), and the (superconducting) order param- 
eter exhibits SO (2) symmetry at low temperature. As 
a result, there is a Berezinskii-Kosterlitz-Thouless phase 
transition to a superconducting state at a finite temper- 
ature rBKTi7^T'?i77 

In the following, we consider only the half-filled case 
where the attractive model maps onto the repulsive 
model studied in the present work. Since the Green's 
function and the spectral function are invariant under 
the particle-hole transformation II79|I , we can directly ap- 
ply the results obtained in the previous sections. The 
phase diagram is shown in Fig. El The crossover lines 
are the same as in Fig. El but their physical meaning 
is different. Below the HF transition temperature , 
the S0(3) order parameter (pq=7r,Aq=o) acquires a fi- 
nite amplitude Aq. This corresponds to the appearance 
of bound particle-hole and particle-particle pairs with a 
size ^0 ~ ^/Aq. Below Tx, directional correlations of the 
order parameter (/Oq=7r, Aq^g) start to grow exponen- 
tially (renormaHzed classical regime) and eventually long- 
range order sets in at the Tc = phase transition. Be- 



cause of the SO (3) symmetry, the ground state can have 
any combination of superconducting and charge-density- 
wave long-range orders. As U increases, the ground state 
smoothly evolve from the BCS to the Bose-Einstein lim- 
its. In the weak-coupling limit {U <C 4t), there is a pseu- 
dogap regime at finite temperature due to the directional 
fiuctuations of the SO (3) order parameter. In the strong- 
coupling limit {U ^ At), between and Tx, there is 
a regime of preformed (local) particle-particle pairs with 
no superfiuid or charge-density-wave short-range order 
(^ ~ 1). Only below Tx do these bosonic pairs begin 
to develop short-range order. At T = 0, the particle- 
particle pairs Bose condense and/or localize, thus giving 
rise to superfiuid and/or charge-density- wave long-range 
orders. The dotted fine in Fig. El is obtained from the 
vanishing of the tunneling DOS p{uj = 0) at zero energy. 



VI. SUMMARY AND CONCLUSION 

We have presented an approach to the 2D half-filled 
Hubbard model which describes both collective spin fiuc- 
tuations and single-particle properties for any value of 
the Coulomb repulsion U. It is valid below a crossover 
temperature Tx where amplitude fiuctuations of the AF 
order parameter are frozen out and AF short-range or- 
der starts to grow exponentially (renormalized classical 
regime) . 

The magnetic phase diagram is obtained from a NLctM 
that is derived from the Hubbard model. The parameters 
of the NLcrM, the bare spin stiffness p° and the spin- wave 
velocity c, are expressed in terms of the mean- value of the 
kinetic energy and current-current correlation functions 
in the HF state. The model is solved by a saddle-point 
approximation within the CP^ representation where the 
Neel field is represented by two Schwinger bosons. Bose- 
Einstein condensation of the Schwinger bosons at zero 
temperature signals the appearance of AF long-range or- 
der. At finite temperature (below Tx), the system is in 
a renormalized classical regime where the AF correlation 
length ^ is exponentially large. The single-particle prop- 
erties are obtained by writing the fermion field in terms 
of a Schwinger boson and a pseudo-fermion whose spin is 
quantized along the (fiuctuating) Neel field. This decom- 
position allows us to approximate the fermion Green's 
function by the product (in real space) of the Schwinger 
boson propagator (which is obtained from the NLcrM) 
and the HF fermionic propagator. 

Our results are summarized in Fig.El which shows the 
phase diagram of the 2D half-filled Hubbard model, and 
Figs. 141101 At weak coupling and zero temperature, our 
theory clearly describes a Slater antiferromagnet with 
an exponentially small AF gap, well-defined Bogoliubov 
QP's carrying most of the spectral weight, and an inco- 
herent excitation background at higher energy. As U in- 
creases, the Slater antiferromagnet progressively evolves 
into a Mott-Heisenberg antiferromagnet with an AF gap 
of order U and a significant fraction of spectral weight 
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transferred from the Bogoliubov QP's to the incoherent 
excitation background. At finite temperature, the Bo- 
goliubov QP's disappear and only incoherent excitations 
survive. Nevertheless, precursors of the zero-temperature 
Bogoliubov QP's show up as sharp peaks in the fermion 
spectral function, with a width of order T. The pres- 
ence of thermal spin fluctuations gives rise to fermionic 
states below the zero-temperature AF gap. At weak cou- 
pling, the latter is completely filled and replaced by a 
pseudogap. The DOS p{'jj) remains however exponen- 
tially small at low energy. At strong coupling and fi- 
nite temperature (0 < T < Tx ~ J), the system is a 
paramagnetic Mott-Hubbard insulator in a renormalized 
classical regime of spin fiuctuations. At higher temper- 
ature, Tx ^ J ^ T < T^^, the system is character- 
ized by the presence of preformed local moments with- 
out AF short-range order. Thus our theory predicts a 
metal-insulator transition at finite temperature between 
a pseudogap phase at weak coupling and a Mott-Hubbard 
insulator at strong coupling. For the 3D Hubbard model, 
we expect a similar phase diagram, but with Tx replaced 
by a true transition line T^ between a paramagnetic phase 
and an AF phase. The weak coupling pseudogap phase 
therefore appears as a consequence of the low dimension- 
ality of the system and the high symmetry (i.e. S0(3)) 
of the AF order parameter. 

At half-filling the attractive and repulsive Hubbard 
models can be mapped onto one another by a canoni- 
cal transformation so that our results also apply to the 
attractive case. AF fiuctuations in the repulsive model 
correspond to q = tt charge and q = pairing fiuctua- 
tions in the attractive model. The corresponding phase 
diagram is discussed in Sec. (see Fig. [T2|l . 

Besides its validity both at weak and strong cou- 
pling, our approach differs from previous weak-coupling 
theorieo-^^.-^^.^^-^'^.^7,-38.39.40.4i,42,4.3,44,4.5,4fi,47 ^^e pseu- 

dogap phase in two respects. First, it takes spin fiuctu- 
ations into account within a highly non-Gaussian theory 
(the NLcrM) and is valid at low temperature (0 < T ^ 
Tx). On the contrary, most of the other approaches as- 
sume Gaussian spin fiuctuations so that their range of 
validity is restricted to T ~ Tx- Second, our NLtrM ap- 
proach is an expansion about the AF ordered state which 
is a valid starting point in presence of AF short-range or- 
der. When calculating fermion propagators, we have to 
consider HF pseudo-fermions interacting with Schwinger 
bosons whose dynamics is determined by the NLtiM. 
Since the HF pseudo-fermions are gapped, we expect a 
perturbative expansion in the pseudo-fermion-boson in- 
teraction to be well-behaved. Our results were obtained 
to lowest order where the fermion Green's function is 
given by the product (in real space) of the HF fermionic 
propagator and the Schwinger boson propagator (which 
is obtained from the NLtrM). This should be contrasted 
with perturbative treatments applied to free fermions in- 
teracting with soft collective fiuctuations where no small 
expansion parameter is available. 

Our NLctM approach is reminiscent of slave-fermion 



theories^SiSiiSS where the fermion is written as the prod- 
uct of a spinless pseudo-fermion and a Schwinger boson 
carrying the spin degrees of freedom. Slave-fermion the- 
ories apply to the t-J model where the Hilbert space is 
truncated by forbidding double occupancy of the lattice 
sites. In our work, the pseudo-fermion also carries a spin, 
which is a necessary condition to describe both the weak 
and strong coupling regimes. 

Our approach bears also some analogies with the work 
of Gusynin et al^^t^^^ on 2D fermion systems with an 
attractive interaction. These authors use a "modulus- 
phase" representation for the SO (2) superconducting 
order parameter which is analog to our "amplitude- 
direction" representation of the SO (3) AF order pa- 
rameter. At low temperature, the phase of the su- 
perconducting order parameter is governed by a SO (2) 
sigma model. The fermion Green's function is calcu- 
lated both above and below the Berezinskii-Kosterlitz- 
Thouless phase transition Tbkt by writing the fermion 
field as the product of a pseudo-fermion and a bosonic 
field which is related to the phase of the order pa- 
rameter. As in our work, a simple decoupling proce- 
dure between pseudo-fermions and bosons is used. A 
pseudogap phase is found both above and below Tbkt- 
Gusynin et al. also point out the necessity to per- 
form a fully quantum-mechanical calculation to describe 
the pseudogap phaseiSi The main difference with our 
work comes from the SO (2) symmetry of the order pa- 
rameter which leads to a finite-temperature Berezinskii- 
Kosterlitz-Thouless phase transition. 

Let us now mention some limitations of our ap- 
proach, (i) The feedback of spin fiuctuations on (pseudo)- 
fermions is not fully taken into account. As a result, we 
miss important effects, like the renormalization of the 
zero-temperature HF gap Aq by quantum spin fiuctua- 
tions. (ii) The crossover temperature Tx, which is iden- 
tified to the HF transition temperature T^^ at weak 
coupling, is overestimated. Due to Kanamori screen- 
ing effects, Tx should be smaller than (iii) 
The NLcrM approach is restricted to low temperature 
(T <C Tx). In particular, it does not give access to the 
crossover regime between the Fermi liquid and the pseu- 
dogap phase at weak coupling. This regime is character- 
ized, as the temperature decreases, by the suppression 
of Landau's QP's. (iv) At finite temperature, we pre- 
dict a metal-insulator transition between a pseudogap 
phase and a Mott-Hubbard insulator. However, being 
a low-energy theory, the NLcrM does not allow to study 
the finite-temperature metal-insulator transition in detail 
(see SecHg. 

But the main shortcoming of our approach is that 
it does not distinguish between Bogoliubov and Mott- 
Hubbard bands. We find a single energy scale (Aq) 
in the density of states p(uj) and the spectral function 
.4(k, w). On physical grounds, we expect instead two 
energy scales, namely Aq and f//2, corresponding to 
Bogoliubov bands (or precursors thereof at finite tem- 
perature) and Mott-Hubbard bands, respectivelyjiSi In 
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the weak coupling limit, Aq depends crucially on the 
nesting properties of the Fermi surface (Slater antifer- 
romagnetism) . On the other hand, the energy scale U/2 
has a purely local origin, which is independent of the 
Fermi surface geometry, and is associated with the Mott- 
Hubbard localization. A proper description of the Mott- 
Hubbard localization would require to treat the charge 
fluctuations beyond the HF approximation for the Ac 
fleld (Sec.nj. In the strong-coupling limit, charge fluc- 
tuations are frozen out. This is the reason why the HF 
saddle point for the amplitude fields Ac and As provides 
an accurate description of the local moments (whose di- 
rection is given by the ilr field) which form in the strong 
coupling limitiSii Note that for U > 4t, Ao U/2 so 
that the system is characterized by a single energy scale. 
At intermediate coupling {U ^ 8t), a four-peak struc- 
ture corresponding to the simultaneous presence of Bo- 
goliubov and Mott-Hubbard bands has been observed in 
numerical simulation a^^i^^ and analytical studiea^^*^ of 
the Hubbard model. Although it misses some aspects of 
the Mott-Hubbard localization, in particular at interme- 
diate coupling, we believe that our theory captures the 
main features of the physics of the 2D half-filled Hubbard 
model. 

There are several directions in which this work could be 
further developed. The most obvious one is to consider 
situations where antiferromagnetism is frustrated due to 
either a non-bipartite lattice or a finite next-neighbor 
hopping amplitude. Doping would also induce magnetic 
frustration. This opens up the possibility to stabilize 
more exotic magnetic orders (e.g. a non-collinear order), 
and/or to reach the quantum disordered and quantum 
critical regimes of the NLcrM (Fig.[2l and study the cor- 



responding fermion spectral functions. 
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APPENDIX A: HF CURRENT-CURRENT 
CORRELATION FUNCTION 

In this appendix we calculate the static uniform 
current-current correlation function 

n;i;:U(j;:(o,o)j;;(o,o))^^. (ai) 

From the definition of the current [Eqs. II25I26II ]. we 
see that its zero-frequency zero-momentum Fourier trans- 
form involved in Eq. ijAlll is given by 

^'^^^''^^ " 7^E^M(k)<i>L'^^<i'k.., (A2) 

where 

vo(k) = 1, (A3) 

v^(k) = 2tsin(fcp), fi^x,y. (A4) 

Using the Wick's theorem to evaluate HF averages of <i> 
fields, we can express ^'f^'^^^, as 



^ z;^(k)v(k')Tr k (^k.'&L.' ) ('J^k'^'^L 



(A5) 



where Tr denotes the trace with respect to the spin indices. Writing the HF propagator [Eq. lfTT|l ] as 

^^k^^-k'^') = (5^,^' [4,k'G(k,Cj) +(5k,k'+7r'73F(k,w)] , 

^iuj — Ck 



G'(k,w) 



and using TT{ai,cr^i) — 2S,y^^>, TT{(Tj,av<yd,cFv') = 2(5^,i/' (2(5^3 — 1) and F(k + 7r,a;) = F(k, cj), we obtain 

' (3N 



i^M^vlk') K,k'G(k,c^)2+,5k,k'+,r(2<5,,3-l)F(k,c.)2] 



k.k' 



(A6) 
(A7) 

(A8) 



nj^^, is thus diagonal in v and v' . One can show that it is also diagonal in ^ and /i'. Indeed, whenever these two 
indices are different, the rhs of ijASp is odd in kx or ky and vanishes after wave-vector summation. Furthermore, 
tio(k + 7r) = wo(k) and v^{\l + tt) = — t;^(k) for fi = x,y, so that 



(3N 



^ ^ «^(k)2 [G(k, ujf + (2<5^,o - 1)(2J,,3 - l)F(k, uof] 
I 



(A9) 



For T <C Tn^, one can perform the Matsubara frequency 
summation in the zero-temperature limit. This gives 

-i5:G(k,.)^ = ^E^(k^-)' 



(AlO) 
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The only non- vanishing correlator functions are therefore 



yrll _ tt22 



n 



n 



00 

33 

yy 



sin^ kx 



(All) 



(A12) 



APPENDIX B: SU(2) GAUGE FIELD 



In this appendix we give a proof of Eqs. (|35l3fi|l . relat- 
ing the Neel and canting fields rir and Lr to the gauge 
field A'j^^ and the rotated canting field lr. Let us recall 
the definition of the gauge field: 



(Bl) 



The index t stands for real-time derivation. Imaginary- 
time results are obtained using r = it. The SU(2)/U(1) 
rotation matrix Rr is defined, up to a U(l) gauge trans- 
formation i?r — > -RrG*""''^^ , by 



cr • rir = RrCT^Rl, 



(B2) 



= i(nrA9^nr)^ 

= \{d,n^?- (B6) 
Using dt = idr we obtain Eq. Ij35|l . Second, recalling that 
the rotated canting vector lr = TZ^^Lr has no component 
along Uz, we can write 



1 



y=l,2 



(rir A a^rir) • Lr, 



(B7) 



hence Eq. l|36|l . 

We now give a derivation of Eq. ljB5|l . The first step 
is to differentiate Eq. IIB2|l . Derivatives of the rotation 
matrix are calculated using Eq. (IB1|I and the identity 
d^Rl = —Rl{df^Rr)Rl which results from the unitarity 
of Rr. We obtain 

(T • d^Ur = -iRrAfj_r<T3Rt + «i?rCr3 Atr-R^ 

= -iRrA^rRlRrCraRl + iRrCrsRlRrAf^rRl- 

(B8) 



which means that the SO (3) element TZr associated to 
Rr maps onto rir. The gauge field A'j^^ is a zero-trace 
Hermitian matrix which can be decomposed on Pauli ma- 
trices cTj,: 



Let us define a new field 



Afj_r — ^ ^'ilr'^f 



u=l,2,3 



(B3) 



where the bold notation denotes the three-component 
vector (Ai,A2,A3). 

The main result of this appendix is the following gen- 



eral form for the A^r field: 



A — 7?^^ 



-rir A a^rir 



-U^ A7^r ^(a^Hr 



(B4) 
(B5) 



Kf^r is some function of position and time, fixed by the 
choice of a gauge. Notice however, that it cannot be any 
function, since it appears in the expression of the gauge- 
field density tensor, which must be zero. 

Eqs. H35I36|I follow quite easily. First, we have 



i'=l,2 



ArJ = i [7^-l(nr Aa^rir)!' 



Afj,r — RrAfj^rR]- 
= Anr ■ a-. 



= 2(7 ■ i A^r A rir 



We have used the identity 

[cr • u, (T • v] 2i(T ■ (u A v). 



(B9) 



Using Afj^r and Eq. ljB2|l we can rewrite Eq. ijBSII as 

cr • a^rir = -i [Af,r, cr • Hr] 

= -i[cr ■ Apr,cr • rir] 



(BIO) 



(Bll) 



Identifying the coefficients of cr in Eq. ljBir)|l and vector- 
multiplying by rir we arrive at 

A^r = ^nr A a^rir + (rir ■ A^r)nr. (B12) 

To conclude, it is sufficient to define the last term in 
Eq. ljB12|l as K^^r and to remark that, owing to the defi- 
nition of A^r, we have A^r = Ti-r-^^r- 
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